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We study the folding of the regular two-dimensional triangular lattice embedded in 
the regular three-dimensional Face Centered Cubic lattice, in the presence of quenched 
random spontaneous curvature. We consider two types of quenched randomness: (1) a 
"physical" randomness arising from a prior random folding of the lattice, creating a pref- 
ered spontaneous curvature on the bonds; (2) a simple randomness where the spontaneous 
curvature is chosen at random independently on each bond. We study the folding transi- 
tions of the two models within the hexagon approximation of the Cluster Variation Method. 
Depending on the type of randomness, the system shows different behaviors. We finally 
discuss a Hopfield-like model as an extension of the physical randomness problem to ac- 
count for the case where several different configurations are stored in the prior pre-folding 
process. 
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The statistical properties of polymerized membranes have been extensively studied 
in the past few years |l],||. Some particular attention was paid to the role of quenched 
disorder in the elasticity of the membrane, with mainly two motivations. The first one is 
to understand the mechanism of the "wrinkling" transition of partially polymerized lipid 
vesicles || . Such membranes undergo a reversible phase transition from a high-temperature 
soft phase with strong fluctuations to a low-temperature rigid and highly wrinkled phase. 
A second motivation is, at a macroscopic level, the study of the statistical properties 
of randomly crumpled paper or, more generally, randomly crumpled elastic sheets 0,||. 
In a random crumpling process, creases are created, which generate random spontaneous 
curvature. When iterated, the random crumpling processes can moreover cause frustration, 
and the crumpled paper may then have many equally probable configurations of minimal 
energy, a usual characteristics of random spin systems. 

Here we study a simple system of two-dimensional polymerized object with quenched 
disorder. As a toy model, we consider the problem of folding of the regular two-dimensional 
triangular lattice in the presence of random spontaneous curvature. Models of folding have 
been introduced in || and studied in |7|||. Originally, the study was restricted to planar 
folding, i.e folding in a two-dimensional embedding space. A more general discrete folding 
model with a three-dimensional embedding space was then introduced and studied in P, 



10.11], describing the folding of the triangular lattice in the regular three-dimensional 



Face Centered Cubic (FCC) lattice. The role of disorder in the planar folding problem was 



analyzed by the authors and Di Francesco in ||12|| .There, disorder was introduced in the 
form of a quenched random bending rigidity. Here we would like to complete our study by 
considering the folding in the three-dimensional FCC lattice in the presence of quenched 
random spontaneous curvature. 



Before we proceed to our study, let us recall the results of [T^] for the planar folding 
problem with random bending rigidity. There, we were interested in modeling the random- 
ness arising from a prior irreversible "crumpling" of the lattice. Such a crumpling results 
in a marking of the lattice bonds with quenched creases on which folds are favored. The 
system can then be described by a Mattis-like model [Ol with Hamiltonian 



T^Mattis = ~K 2J TiTjOiOj , (1) 

n.n.(ij) 

where the variables o~i = ±1 describe the (up or down) normal to the triangle i in the 
folded configuration and n.n.(ij) means summation over nearest neighbour pairs. The 
disorder variables r% = ±1, accounting for the prior irreversible marking, define a random 
bending rigidity Kij = KtiTj, and are "frozen" according to a specified probability dis- 
tribution |12| . For the variables a to represent actual folded configurations of the lattice, 



the six neighboring spins on an elementary hexagon of the lattice, o~i{i = 1, 2, • • • , 6), must 
satisfy the "physical" constraints ||||: J2i=i°~i = mod 3. This condition in particular 
prevents from absorbing the disorder in a simple change of o~i into a^. Similarly, if the 
corresponding disorder variables Tj(z = 1, 2, • • • , 6) arise from a pre-folding process, they 
should also obey the physical rule Yui=i r « = m °d 3. This condition on the r-variables 
is essential for the system to develop a large K "frozen phase" where the prior irreversible 
folded shape is recovered. If the r variables are free ±1 variables which do not satisfy 
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the physical constraint above, the system becomes frustrated and the lattice remains in a 
disordered phase. 

We now would like to extend these results to the case of the three-dimensional FCC 
folding problem. A folding of the triangular lattice in the FCC lattice is simply a mapping 
sending each vertex of the triangular lattice onto a vertex of the FCC lattice, with the 
requirement that neighboring vertices on the triangular lattice remain nearest neighbors 
in the FCC lattice [^|,|T0|] , i.e. belong to the same triangular face. The FCC lattice is 
indeed made of octahedra and tetrahedra in contact by their triangular faces. Elementary 
triangles of the triangular lattice are thus sent onto elementary triangular faces of the FCC 
lattice. In the folded configuration, two adjacent triangles can form some relative angle 9, 
with one of the four following values: 
(i) 9 = 180° — no fold: the triangles are side by side; 
(ii) 9 = 0° — complete fold: the triangles are on top of each other; 
(iii) 9 = arccos(l/3) ~ 71° -- fold with acute angle: the two triangles lie on two adjacent 

faces of the same tetrahedron in the FCC lattice, and 
(iv) 9 = arccos(— 1/3) ~ 109° -- fold with obtuse angle: the triangles lie on two adjacent 

faces of the same octahedron in the FCC lattice. 
It was shown in [|| that these four types of folds can be understood as the superposition 
of the domain walls of two Z2 variables a = ±1 and z = ±1 living on the faces of the 
triangular lattice. The relative values Act = o"20"i and Az = z±Z2 for two neighboring 
triangles indicate which type of fold they form, with the correspondence displayed in 
Table I. 



Act 


Az 


9 


angle 


1 


1 


180° 


no fold 


-1 


1 


0° 


complete fold 


1 


-1 


71° 


acute fold 


-1 


-1 


109° 


obtuse fold 



Table I: The relative folding state of two neighboring triangles according to 
the relative values Act and Az. 

In order to describe an actual allowed folded state, the a and z variables are subject to 
two basic folding rules involving the values Oi and z% (i — 1, • • • , 6) on the six neighboring 
triangles forming an elementary hexagon in the lattice. For each hexagon, the variables a 
must satisfy the following first folding rule: 



y Oi = mod 3. 



(2) 



i=i 



This rule is identical to the rule of planar folding [[7] although its interpretation here is 
slightly different 0,^O|. A second basic folding rule involves both the z and a variables 



and reads: 

i 

I 1 ZiZ i+ i = 1 for c = 0, 1, 2 ; 1(c) = {i : 2_. °k = c mod 3}. (3) 

iei(c) fc=i 

With the "physical" constraints (|j) and (0) , one finds exactly 96 possible hexagonal config- 
urations for the six triangles surrounding any of the vertices of the triangular lattice. Note 
that the planar folding problem can be recovered by freezing the z variable to z% = +1 
globally for all triangles. One is then left with exactly 11 possible hexagonal configurations. 
In the absence of disorder, the folding energy is simply -Epure = — K cos (9) per lattice 
bond, with K the bending rigidity parameter. In terms of the variables CTi, a 2 and z±, z 2 
of the two triangles forming the fold, the folding energy simply reads: 

K 

E vure = -—a 1 a 2 (l + 2z 1 z 2 ). (4) 

The total folding energy is the sum of all elementary folding energies for all the bonds of 
the triangular lattice. 

Disorder can be put in the model by introducing quenched disorder variables Tj and Wi 
describing the pre-folded state created by the irreversible crumpling process. The domain 
walls for the variables r and w encode the four possible types of created creases with angle 
ip = 180°, 0°, 71° or 109° according to a table similar to Table I. Of course, in order to 
describe an actual pre-folded state, the r and w variables are subject to two local folding 
rules similar to (|2]) and (|3|). The presence of random creases directly leads to a random 
spontaneous curvature in the system, encoded in the angle ifi of the crease. Given this 
angle, the energy becomes minimum when the angle 9 of the fold is such that 9 — ip. We 
shall thus consider the following bending energy: 

E--K cos(9-^), (5) 

where K measures the strength of this bending energy and where the quenched random 
variable ip describes the quenched random spontaneous curvature in the system. As men- 
tioned above, the variables 9 and ij) take four values, leading to 16 possible fold/crease 
configurations. We would like to express the energy (|5|) in terms of the spin variables Oi, 
Z{, Ti and Wi (i — 1,2) on the two neighboring triangles on each side of the fold, as we did 
in (U) for the pure case without disorder. The relative values Act = o"iO"2 and Az = z\z 2 
on one hand and those of At = t\t 2 and Aw = w±w 2 on the other hand fix the angles 9 
and ij) and thus the energy (|5|). 

We can make use of the symmetry 9 <-» ijj in (|5|) (i.e. (Act, Az) ^ (At, Aw)) and of 
the symmetry (9, ip) <-> (180° - 9, 180° - ip) (i.e. (Act, At) <-> -(Act, At)) to ensure that 
the bending energy has the two independent symmetries At <^ Act and Aw ^ Az, and is 
even in (Act + At)/2, leading to the general form (for Z 2 variables): 

E (Az + Aw) , .. . . , ,. . N 

— . = e + a- + 6i(ActAt) + b 2 (AzAw) 

2 fr.\ 

,. . .(Az + Aw) ,, . . w . . . 
+ c(ActAt)^ + ^(ActAt)(A^Aw) 



Fold Config. 


Disorder Config. 


Bending Energy/X 


No fold 


No Crease 


-1 


No fold 


Acute Crease 


1/3 


No fold 


Complete Crease 


1 


No fold 


Obtuse Crease 


-1/3 


Acute Fold 


Acute Crease 


-1 


Acute Fold 


Obtuse Crease 


-7/9 



Table II: The folding and crease state of two neighboring triangles and the 
corresponding bending energy /K for the 6 independent configurations. 

involving 6 constant coefficients e, a, bi, 62, c, d to be determined hereafter. These coeffi- 
cients are simply obtained from the values of the 6 independent (i.e. not related by the 
above symmetries) pairs of folding and disorder configurations given in table II. 
From this table, we determine the bending energy to be: 



E = — 



K 
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_ 2 + 4 ( Az + Aw ) _ ( AaAr ) _ 2(AzAw) 



4(AaAr) ^ Az + A ^ - 4(AoAt)(AzAw) 



(7) 



K 



[Aa(l + 2A;z)Ar(l + 2Aw) + 2(1 - Az)(l - Aw)} 



The total bending energy is again the sum of all elementary folding energies for all links 
of the triangular lattice. It is interesting to check several limiting cases of the above 
formula. The pure folding problem, without disorder, can be recovered by constraining 
the disorder variables according to At = Aw = 1. The bending energy (|7|) reduces to 
E = —KAa(l + 2Az)/3, i.e to Eq.(Q). The case of planar folding with disorder is obtained 
by setting Az = Aw = 1 in (^) and the bending energy reduces then to E = —KAtAo, 
as in Eq. (|T|) . 

In addition to the bending energy (0), we also introduce an external field H r associated 
with the variable (ar+wz), which is a rough measure of how close the folding configuration 
(a, z) is from the disorder configuration (r, w). This definition of the external field is not 
canonical and many other definitions are equally acceptable. The main motivation for 
introducing this external field is technical, i.e. the necessity to first prepare a solution 
with explicit broken symmetry to be able to eventually reach a solution with spontaneous 
symmetry breaking by tuning H r to zero. Our total Hamiltonian is thus given by: 



n 



K 



— ^J [tjTj (l + 2wiWj)<Ji(jj (l + 2zjZj) + 2(1 -WiWj)(l-z i Zj)]-H r y^ j (cTiTi + 



ZiWi 



n.n.(ij) 



(8) 



To analyze the properties of the system, we will consider different order parameters. In 
view of what we know about the pure system without disorder, it is useful to divide the 
original triangular lattice into the two subsets A and B made of all triangles pointing up 
and all triangles pointing down respectively in the original fiat triangular lattice. This in 
particular allows us to define "staggered" average values. We will be interested primarily 
in the four following average values: 



Sa e 


= <">A 


Z A = 


= <z)a 


F x = 


■ (or) 



(9) 

F 2 = (zw) , 

where the brackets denote the average over the configurations at fixed disorder for an 
arbitrary given triangle (taken moreover in the subset A when we add the index A to the 
brackets) and the over-line denotes the average over the quenched disorder. One has of 
course — 1 < Sa, Za, Fi, Fi < 1. Non-zero values of the "frozen" order parameters F\ 
and F2 indicate that the membrane is trapped in the configuration given by the disorder 
variables r and w. The reason why we restrict ourselves to the subset A in Sa and Za 
is that, for the "pure" system, that is the model without disorder, and at K = 0, it was 
established that the lattice is found in a phase where Sa — —Sb ~ 0.874560 and Za = 
Zb = 0. This can be interpreted as a strong preference for the lattice to wrap on octahedra 
in the FCC lattice. In this phase, clearly the full average value S = (Sa + £b)/2 = is 
not a good quantity and the correct order parameter is S st = (Sa — S-q)/2, or more simply 
Sa itself. 

We now come to the question of the precise form of the probability distribution for 
the quenched disorder variables r and w. As has been discussed previously, these disorder 
variables should obey the two folding constraints in order to correspond to some prior 



folding of the lattice. As in our previous work in [pL2| , we can take advantage of the 
solution of the pure system and simply assume that the disorder distribution is described 
by a particular equilibrium distribution of this pure system. Since there is no physical 
reason to introduce an energy scale in the distribution of the disorder variables and because 
we want to treat as equiprobable all pre-folded configurations, the natural choice is to take 
the distribution of the pure system at K — 0. As we just mentioned, the triangular 
lattice is then in a "octahedrally" folded phase. The disorder configuration will thus also 
have the same nature, i.e. fst = 0.874560. This means that the disorder configuration is 
dominated by obtuse and complete creases. We will refer to the disordered model with this 
particular probability distribution as Model 1. For comparison, we will also study a model 
without the physical constraints on the r and w variables, i.e a case were the two random 
variables take ±1 values with equal probability, independently on each triangle. All the 
(2 x 2) 6 = 4096 hexagonal disorder configurations are then possible and equiprobable. We 
shall call refer to this second model Model 2 and will compare its behavior with that of 
Model 1. 

To analyze the properties of both Models 1 and 2, we rely on the same method we 
used in |12[ for the planar case, i.e. the Cluster Variation Method (CVM) in its hexagonal 



approximation. The CVM uses a variational principle on the free energy of the system 
together with a suitable truncation of the cumulant expansion of the entropy at the level 
of some maximal clusters (here the hexagons) [fn| - |T6|| . It is applicable to the statistics of 
both pure systems without disorder and to systems with quenched random disorder. In 
this case, it allows in particular to evaluate average values such as (H). This method has 
been first used for the planar and FCC folding problems in [17]. Although it is only an 
approximation since the entropy is evaluated from a truncation of its cumulant expansion, 
it has been shown in the most simple cases of folding problems that its results compare 
extremely well with exact predictions ]T7J . For more complex cases where no exact results 
exist, we still believe that this method allows for reliable and accurate predictions, as far 
as the nature of the phases and of the transitions are concerned. We refer to our previous 



work [|T2|] for a detailed description of this method and its implementation for disordered 
systems. 
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Fig.l: The free energy and the order parameters Sa (dashed), Zx (solid), 
Fi } c (dashed) and F 2 (solid) as a function of the bending energy K and the 
external field H r . These results are for Model 1, i.e. with disorder variables 
satisfying physical constrains. 

We now present the results of our CVM analysis. Our main result is, as for the planar 
system, the dependence of the existence of a frozen phase Fj ^ (i = 1, 2) on the type 
of disorder. In the case of Model 1, a discontinuous transition occurs at K = Kpz ~ 0.44 
from a low K octahedrally folded phase to a large K frozen phase. Model 2 does not 
develop such a phase and the lattice is always found in the octahedrally folded state (this 
was checked up to K = 10.0). We have also looked at the effect of the external field H r 
alone. In this case, both systems show a discontinuous transition to a partially frozen 
phase < .Fi < 1, with however very different resulting states. In the frozen phase of 
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the Model 1, the Fi almost saturate to 1 and the lattice is thus completely trapped in 
the configuration specified by the disorder variables. In the case of Model 2, the values of 
Fi are much smaller, which means that the degree of freezing is far from complete. Such 
a different character comes from the frustration appearing in the system for Model 2, a 
result similar to what we found in the planar folding model [12]. The external field tries to 



put the lattice configuration into a disorder configuration which is in general not accessible 
due to the folding constraints on the a and z variables. In other words, there is not a 
unique ground state for arbitrary disorder variables. 
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Fig. 2: The free energy and the order parameters 5a (dashed), Zx (solid), 
Fi } c — Fi (dashed) and F 2 (solid) as a function of the bending energy K 
and the external field H r . These results are for Model 2, i.e. with disorder 
variables completely random. 



In Fig.l and Fig. 2, we show the behaviors of the free energy and the different order 
parameters Sa, Za, F± t c = F± — Sa^a and F 2 for both the Model 1 and the Model 2 as 
a function of the bending energy K and the external field H r . We display F\c instead 
of Fi, because in the case of Model 1, F\ takes non-zero values even at (K, H r ) = (0, 0). 
Indeed, when K = H r = 0, the disorder variables r and w are decoupled from the physical 
degrees of freedom a and z. From the definition of F\, its value is then simply given by 
the product of the value Sa of the pure system at K = and the quantity ta, both equal 
to 0.874560, leading to F\ = (0.874560) 2 ~ 0.765. This non-zero value is subtracted if we 
use the "connected part" F^c instead. 

In the case of Model 1, by increasing the bending rigidity K, the systems undergoes 
a phase transition from the octahedrally folded phase (Sa ^ 0, Za = 0, -Fie 7^ and 
F 2 = 0) to the frozen phase (Sa ~ tX = 0.874560, Z A = 0, F hC ~ l.-(0.874560) 2 ~ 0.235, 
F 2 ~ 1). The values of F{ saturate to 1, meaning that the system is almost completely 



trapped in the configuration specified by the disorder variables. The external field also 
induces a phase transition. In this case, a weak field already creates a non-zero value of 
i*2 and the lattice is thus in a weakly frozen phase. A discontinuous transition occurs to 
a highly frozen phase at with strong freezing Fi ~ 1 (i = 1, 2). 

In the case of Model 2, increasing the bending rigidity does not give rise to a frozen 
phase (in the range K < 10.0 at least were we performed our analysis). The octahedral 
order 5a 7^ persists for all K and the triangular lattice remains in the octahedrally 
folded phase. With the external field H r , the octahedral order disappears abruptly and a 
partially frozen phase Fi 7^ (i = 1, 2) appears at H rj c = 0.18. Still, the values of Fi and 
F2 are small in this phase compared with the large H r frozen phase of Model 1. 



In the present work and in our previous work |[12|| , we have studied triangular lattice 
folding models with random bending rigidity and random spontaneous curvature. These 
models can be seen as toy models for crumpled paper in the sense that they combine the 
effect of disorder and of geometrical metric constraints on a geometrical two-dimensional 
object. We have in mind the situation where the paper is crumpled only once and the 
system thus does not have frustration and can recover this crumpled state. A natural 
question is now what happens if the paper is crumpled several times. In this case, we expect 
that the lattice should store several creases configurations. Let us discuss now how this can 
be implemented in our model. For simplicity we return to the planar folding problem with 
random bending rigidity (see Eq.(|l|)). We can imagine that the lattice, after p crumpling 
processes, stores p different configurations, which we denote by rf (/x = l,---,p). Then 
one possible choice is the following bending energy, written by analogy with the Hopfield 
model 0@: 

^; = f£rfTf. (10) 

We note however two differences with the usual Hopfield model. In our problem, the 
spin variables Oi and the disorder random variables rf correspond to folded configurations 
of the lattice, and must satisfy the planar folding constraint (Eq.([|)). Moreover, in the 
Hopfield model, the interaction is long-ranged {Yla j))i while in our model, it is simply 
short-ranged (E n .n.(y))- 

Disorder 

Configurations K . 






(a) (b) 

Fig. 3: An elementary hexagon with two different memories r 1 and r 2 . We 
show on the r.h.s. the corresponding values of Kij according to (a) Eq.(10) 
and (b) Eq.(ll). 



The choice ([K]) for the bending energy is however somewhat unsatisfactory. To see why, 
we consider an elementary hexagon which tries to store the 2 disorder configurations of 
Fig. 3. With the definition of Eq.([K|), the bending rigidity K{j at the creases becomes zero 
(see Fig.3-a), a result which is somewhat unrealistic. We would instead expect that, if one 
folds the hexagon as in the Fig. 3, the bending rigidity becomes negative at all the creases 
and remains positive where there is no crease (see Fig.3-b). A more natural choice for the 
bending energy is then 

K ij =Kmin[T?T?(n = l,---,p)]. (11) 

This choice corresponds to a complete irreversibility of the process of marking creases. 

Disorder Lowest Energy Configurations 

Configuration 





Fig. 4: An elementary hexagon with unphysical disorder and the correspond- 
ing lowest energy configurations satisfying the folding constraint. 

With such a disorder, we now find for the case of Fig. 3 seven ground states in com- 
petition (Fig. 4). This situation is the exactly similar to the case of unphysical disorder, 
where the r variables do not satisfy the folding constraint. In our previous work, we 
have shown that the strong frustration caused by such an unphysical disorder prevents the 
system from having a frozen phase at large K. Based on this result, we could conclude 
that, in the presence of many different stored configurations, the system remains always 
disordered. On the other hand, as we increase the number of irreversible crumplings, the 
number of creases becomes large and the system becomes almost identical to a pure model 
with negative bending rigidity. Such a pure system is known to be at large K in a so-called 
"piled up" ordered phase. We can assume that one develops such an order if the number 
p of stored configuration exceeds a critical value p c i- Finally, from the properties of the 
Hopfield model, on can also imagine the system to be in the so-called spin-glass phase. 

From the above discussion, we can propose the following conjecture on the statistical 
properties of the above model at K = oo, which is of particular importance because 
crumpled paper can be considered as an infinite elastic constant limit K — > oo (or T = 0) 
of the lattice system. The triangular lattice with p stored configurations could be in the 
retrieval phase for very small p < p c ±, where the system can almost store and retrieve p 
patterns. For p c \ < p < p C 2, the system would be in a spin-glass phase or in a disordered 
phase. Finally, for p > p c2 , the system would be in the piled up phase. The values of 
Pd and p C 2 are expected to be very small, because the interaction is short-ranged and the 
probability that each elementary hexagon has two or three creases is very high even with 
one stored configuration. 

9 



Finally, we also expect that the many crumpled configurations which are caused by 
successive crumplings of the paper are not independent from one another and we can 
suppose that they look quite similar. This weakens the frustration, which might keep the 
paper in the retrieval phase for larger p, i.e. increase the value of p& and p C 2- In fact, it 
was also reported experimentally || that the paper has a good memory of the previous 
crumpled configuration. Even in the retrieval phase, the energy landscape in the phase 
space of the system can be very complex. It is this complex nature which might appear 
as the universal power law of the noise emitted from the crumpled elastic sheets || . The 
thermodynamics of the above model are left for a future study. In addition, more realistic 
bending energy could be considered. To conclude, we can say that the problem of multiple 
stored configurations is still very open. 
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